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y '. Abstract 



We analyze the hadronic freeze-out in ultra-relativistic heavy ion colhsions at 
RHIC in a transport approach which combines hydrodynamics for the early, 
dense, deconfined stage of the reaction with a microscopic non-equilibrium 
model for the later hadronic stage at which the hydrodynamic equilibrium 
assumptions are not valid. With this ansatz we are able to self-consistently 
calculate the freeze-out of the system and determine space-time hypersurfaces 
for individual hadron species. The space-time domains of the freeze-out for 
several hadron species are found to be actually four-dimensional, and differ 
drastically for the individual hadrons species. Freeze-out radii distributions 
are similar in width for most hadron species, even though the il.~ is found to 
be emitted rather close to the phase boundary and shows the smallest freeze- 
out radii and times among all baryon species. The total lifetime of the system 
does not change by more than 10% when going from SPS to RHIC energies. 



Ultra-relativistic heavy ion collisions are the only means available to investigate highly 
excited dense nuclear matter under controlled laboratory conditions. In such collisions it 
is sought to recreate a Quark Gluon Plasma (QGP), the highly excited state of primordial 
matter which is believed to have existed shortly after the creation of the universe in the Big 
Bang (for recent reviews on the QGP, we refer to [||]). 

Transport theory has been among the most successful approaches applied to the theoreti- 
cal investigation of relativistic heavy ion collisions. Microscopic transport models attempt to 
describe the evolution of the heavy ion reaction from some initial state up to the freeze-out of 
the newly produced particles on the basis of elementary interactions. The basic constituents 
in such models are either hadrons [0,0 or partons [Q. At RHIC energies, however, both, 
partonic and hadronic, degrees of freedom might be equally important and have both to 
be treated explicitly ^. However, in such microscopic transport models, the QG-matter 
to hadron matter transition, i.e. the hadronization stage, has to be modeled in an ad-hoc 
fashion, whereas hydrodynamic approaches ||6|-[10| incorporate this as a phase transition. 



This can be done in a consistent way, respecting the laws of thermodynamics (which is not 
always the case in microscopic transport models). The drawback of hydrodynamics, how- 
ever, is that in the later reaction stages the basic hydrodynamical assumptions break down. 
For the freeze-out of the system a decouphng (freeze-out) hyper-surface must be specified 
(or fine-tuned to existing data). 

In this letter, we use boost-invariant hydrodynamics to model a first order phase tran- 
sition from a QGP to a hadronic fluid, and combine it with a non-equilibrium microscopic 
transport calculation for the later, purely hadronic stages of the reaction. With this ansatz 
we are able to self-consistently calculate the freeze-out of the system: no decoupling hyper- 
surface is imposed by hand, but the space-time points are rather determined by an interplay 



between the (local) expansion scalar du [11,10| (where u is the collective flow four- velocity). 



the relevant elementary cross sections, and the equation of state (EoS), which actually 
changes dynamically as more and more hadron species decouple. 

Let us first briefly describe the hydrodynamical model employed here: For a more detailed 
discussion we refer to refs. P, |l2| , p!3| . For simplicity, boost-invariant longitudinal flow is 
assumed. For ultrarelativistic collisions, this should be a reasonable first approximation in 
the central rapidity region. Cylindrically symmetric transverse expansion is superimposed. 
For T > Tq = 160 MeV the well-known MIT bag model equation of state |]14| is used, 
assuming for simplicity an ideal gas of quarks, anti-quarks (with masses ruu = ma = 0, 
ms = 150 MeV), and gluons. For T < Tc an ideal hadron gas is employed that includes 
the complete hadronic spectrum up to a mass of 2 GeV. At T = Tc, (/is = /^s = 0) we 
require that both pressures are equal, which fixes the bag constant to B = 380 MeV/fm^. 
By construction the EoS exhibits a first order phase transition (which is also expected in 
QCD for the quark-hadron phase transition in the case of three quark flavors). 



The model reproduces the measured px- and my-spectra of hadrons at the SPS, when 
assuming that hydrodynamic flow sets in on the proper time hyperbola 11 = 1 fm/c P,p!2 



This is a value conventionally assumed in the literature, cf. e.g. |^. Due to the higher parton 
density at midrapidity, thermalization may be reached earlier at RHIC [^ . As in refs. P,p!2 



we assume here r, = Rt/10 = 0.6 fm. The effects of variations of Tj and Tc will be discussed 
in a future publication |TB|. Moreover, we use the initial average energy and baryon densities 
e(rj) = 20 GeV/fm'^ and p(tj) = 2.3pQ, which lead to dNs/dy = 25 and s/p^ = 205 (a bar 
symbols an average over the transverse plane). The initial energy and net baryon densities 
are assumed to be distributed in the transverse plane according to a so-called "wounded 



nucleon" distribution oc |a/1 — rj^/Rj^, with transverse radius Rt = 6 fm. For this set of 
parameters, the initial transverse energy at midrapidity is dET/dy = 1.3 TeV. Due to the 
work performed by the isentropic expansion, it decreases to 716 GeV on the hadronization 
hypersurface. The microscopic treatment of the hadronic dynamics following hadronization 
(see below) yields dEx/dy = 714 GeV at kinetic freeze-out. Thus, the late hadronic evolution 
at RHIC energy is not isentropic. 

After specifying the initial conditions and the EoS, we determine numerically the hydro- 
dynamical solution between the t = Ti hyperbola and the hadronization hypersurface, where 
we apply the Cooper-Frye formula [|T6| to obtain the hadron spectra. However, in contrast 
to the usual procedure we do not integrate over the hypersurface, because we further on also 
need the space-time distribution of hadrons emerging from the mixed phase, not only their 
momentum-space distributions |[T^. The ensemble of hadrons thus generated is then used 
as initial condition for the non-equilibrium microscopic transport model Ultra-relativistic 
Quantum Molecular Dynamics (UrQMD) ||3|JT7|1. The UrQMD model contains hadronic 
(and string) degrees of freedom - all hadronic states can be produced in string decays, s- 
channel collisions or resonance decays. Tabulated and parameterized experimental cross 
sections are used when available. Resonance absorption, decays and scattering are handled 
via the principle of detailed balance. The UrQMD model has been extensively tested in the 
SIS, AGS and SPS energy domain and provides a robust description of hadronic heavy-ion 
physics phenomenology. An extensive description of the model, as well as comparisons with 
various available data can be found in P,p!7|. 

During the mixed phase the system is either described locally within the hydrodynamical 
framework (as long as a non-zero fraction of the fluid in the cell consists of quark and gluons) 
or within the microscopic transport (in the case of pure hadronic matter). Therefore there 
exists a time interval during the reaction in which both models are applied in parallel, even 
though they never refer locally to the same space-time volume. 

Let us now turn to the reaction dynamics of central (impact parameter 6 = fm) Au+Au 
collisions at RHIC energies {^/s = 200 GeV per incident coUiding nucleon-pair). We start 
with the analysis of the freeze-out hypersurfaces of pions and nucleons, the most abundant 



meson and baryon species in the system, restricting ourselves to the central rapidity region 
y = ll/CA/l < 0.5. Figure shows the freeze-out time distributions and the transverse 
radius distributions for both, pions and nucleons. The top row shows the result of a pure 
hydrodynamical calculation up to complete hadronization, with subsequent hadronic decays, 
but without hadronic reinteraction. The bottom row shows the same calculation with full 
microscopic hadronic dynamics added. 

The freeze-out characteristics of both, pions and especially nucleons, are significantly 
modified due to the hadronic interaction phase. The average transverse freeze-out radius 
of the pions changes from 7.8 to 9.5 fm and that of the protons doubles from 5.4 to 11.3 
fm. Their respective average freeze-out times change from 17.2 to 23.1 fm/c (pions) and 
from 11.3 to 25.8 fm/c (protons). As the meson multiplicity in the system is fifty times 
larger than the baryon multiplicity, baryons propagate through the relativistic meson gas 
- they may act as probes of this highly excited meson medium. Thus, a first estimate of 
the duration of the hadronic phase is Ar ^ 13 fm/c. Its transverse spatial extent is on the 
order of Ar^ ~ 6 fm. 

The Hydro+UrQMD model predicts a space-time freeze-out picture which is drastically 
different from that usually employed in the hydrodynamical model, e.g. in refs. i8 HTll , [T9| , ^ : 
Freeze-out here is found to occur in a four- dimensional region within the forward light- 
cone |21[] rather than on a three-dimensional "hypersurface" in space-time [|I6|. Similar 
results have also been obtained within other microscopic transport models [T^ when the 
initial state was not a quark-gluon plasma. This finding seems to be a generic feature of 
such models: the elementary binary hadron-hadron interactions smear out the sharp signals 
to be expected from simple hydro. This predicted additional fourth dimension of the freeze- 
out domain could affect the HBT parameters considerably. 

This does not mean that the momentum- distributions alone can not be calculated assum- 
ing freeze-out on some effective three-dimensional hypersurface. (For example, if interactions 
on the outer side of that hypersurface are very "soft" , the single-particle momentum dis- 
tributions will not change anymore, while the two-particle correlator does change. Thus, 
the freeze-out condition, e.g. the temperature, as measured by single-particle spectra and 
two-particle correlations [^ needs not be the same.) 

The shapes of the freeze-out hypersurfaces (FOHS) show broad radial maxima for in- 
termediate freeze-out times. Thus, transverse expansion has not developed scaling-flow (in 
that case the FOHS would be hyperbolas in the t — rT plane). Moreover, the hypersur- 
faces of pions and nucleons, and their shapes, are distinct from each other (as also found in 



i3| JT0|JT8|j23[| at the lower BNL-AGS and CERN-SPS energies). Thus, our calculation contra- 



dicts the ansatz of a unique freeze-out hypersurface for all hadrons, cf. also refs. |T^,[T^ 

Figure ^ shows the transverse freeze-out radius distributions for vr, K, p, A, S and Q~ . 
They are rather broad and similar to each other, even though the Q~ shows a somewhat 



narrower freeze-out distribution. The average transverse freeze-out radii are 9.5 fm for pions, 
10.2 fm for kaons, 11.3 fm for protons, 11.6 fm for Lambda- and Sigma-Hyperons, 14.2 fm 
for Cascades, but only 7.3 fm for the Q~ . The freeze-out of the Q~ occurs rather close to the 
phase-boundary |12[ , due to its very small hadronic interaction cross section. This behavior 



could be responsible for the experimentally observed hadron-mass dependence of the inverse 



slopes of the mr-spectra at SPS energies ||2^. For the Q~ , the inverse slope remains practi- 
cally unaffected by the purely hadronic stage of the reaction, due to its small interaction cross 
section, while the flow of p's and A's increases [l^. By comparing the transverse freeze-out 
radii of the hydrodynamical calculation (up to hadronization, including subsequent hadronic 
decays, but no hadronic reinteractions) with the Hydro-|-UrQMD calculation, which include 
microscopic hadronic dynamics, the thickness Avhad of the hadronic phase can be estimated 
by computing the difference: Arhad = {ftjr^"^ ) ~ (''"tj/"^ ^ ecays^^ These values for 

Avhad are: 1.7 fm for pions, 3.1 fm for kaons, 5.8 fm for protons, lambda- and sigma-hyperons 
as well as cascades and 2.6 fm for the Q~ . 

Another issue of interest is the predicted significant increase of the lifetime of the system 
from SPS to RHIC energies |^. Figure ^ shows that in our model, which exhibits a first 
order phase transition, there is between SPS and RHIC no difference in the freeze-out time 
distributions of vr, p, and Q^\ Origin of this prediction is that we include many more states in 



the hadronic EoS, which speeds up hadronization considerably P,pO|. Furthermore, decays 



of resonances (which were not treated in P|) mask the remaining small increase of the 
hadronization time. Note that the multistrange Q~ baryons freeze out far earlier than all 
other baryons, as discussed already previously in the context of figure ^ The duration of 

4-1, V, J • • -i- 4-- t, A / Hydro+UrQMD\ / Hydro+had. decays \ • i 

the hadronic remteraction phase, AThad = Vfr ) ~ V/r / remams nearly 

unchanged, e.g. at 5.9 fm/c for pions, 8.0 fm/c for kaons, 14.5 fm/c for protons, 15.4 fm/c 
for hyperons and 8.0 fm/c for the Q~ . 

Note that the lifetime of the prehadronic stage in this approach is a factor of 2 — 3 
longer than when employing the parton cascade model (PCM) 0,^ for the initial reaction 
stage. It will be interesting to check whether this is related to the first-order phase transition 
built into the EoS which is used here. The final transverse freeze-out radii and times (after 
hadronic rescattering) , however, are very similar in both approaches ||^. 

So far, we have only discussed the kinetic freeze-out of individual hadron species, which 
is the most precisely determinable freeze-out quantity of the system. However, apart from 
the kinetic freeze-out, the chemical freeze-out of the system, which fixes the chemical com- 
position is of interest. The top frame of figure ^ shows the time evolution of on-shell hadron 
multiplicities. The dark grey shaded area indicates the duration of the QGP phase whereas 
the light grey shaded area depicts the mixed phase (both averaged over r-^; only hadrons 
that have already "escaped" from the mixed phase are shown). Hadronic resonances are 
formed and are populated for a long time (~ 20fm/c). When the mixed phase ceases to 



exist, the hadron yields have not yet saturated (even if resonance decays are taken into 
account). This is due to inelastic hadron-hadron collisions. In particular the yield of anti- 
protons drops strongly - more than 60% of the baryon-anti-baryon annihilations occur after 
the phase-coexistence is over (c.f. the lower frame of this figure). The yields of all stable 
hadrons saturate at approximately 25 fm/c. Only then may the system be viewed as chem- 
ically frozen-out. Since resonance decays have not been included into our estimate of the 
saturation time, this number may be viewed as an upper estimate of the chemical freeze-out 
time. 

By comparing different final hadron yields resulting from the hydrodynamical calculation 
(up to hadronization, including subsequent hadronic decays, but no hadronic reinteractions) 
to that of the Hydro-|-UrQMD calculation, which includes microscopic hadronic dynamics, 
we can quantify the changes of the hadrochemical content due to hadronic rescattering: 
especially the multiplicities of (anti)baryons vary at least by 10%, those of protons and 
antiprotons even up to 30% (vr : +9.3%, K : -5%, Y : +12%, p : -21%, p : -31%, Y : 
+11%). These changes clearly indicate that in our model chemical freeze-out of (anti- 
)baryons and (anti-)hyperons does not occur directly at the phase-boundary. Also, note 
that (unlike in ideal chemical equilibrium) baryon number is "shuffled" from non-strange to 
strange baryons. 

The bottom frame of figure ^ shows the rates for hadron-hadron collisions. Meson- 
meson (MM) and - to a lesser extent - meson-baryon (MB) interactions dominate the 
dynamics in the hadronic phase. However, the BB collisions outnumber BB reactions, in 
clear contrast to SPS. This is a consequence of the fact that the BB annihilation cross 
sections at small relative momenta increase faster then the total BB cross section [^ . In the 
case of (approximate) baryon-antibaryon symmetry, one therefore expects more BB than 
BB interactions, as seen in figure |. 

All collision rates reach their maxima at the end of the mixed phase - then they decrease 
roughly according to a power-law. After k, 35 fm/c, less than one hadron-hadron collision 
occurs per unit of time and rapidity - at this stage the system can be considered as kinetically 
frozen-out. 

In summary, we have analyzed the hadronic freeze-out in ultra-relativistic heavy ion 
collisions at RHIC in a transport approach which combines hydrodynamics for the early, 
dense, deconfined stage of the reaction with a microscopic non-equilibrium model for the later 
hadronic stage at which the hydrodynamic equilibrium assumptions are not valid anymore. 
Within this approach we have self-consistently calculated the freeze-out of the hadronic 
system and accounted for the collective fiow on the hadronization hypersurface generated 
by the QGP expansion. 

We find that the space-time domains of the freeze-out for the investigated hadron species 
are actually four-dimensional, and differ drastically between the individual hadrons species. 
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The thickness of the hadronic phase is found to be between 2 fm and 6 fm, depending 
on the respective hadron species. Its lifetime is between 5 fm/c and 13 fm/c, respectively. 
Freeze-out radii distributions are similar in width for most hadron species, even though 
the Q~ is found to be emitted rather close to the phase boundary and shows the smallest 
freeze-out radii and times among all baryon species. The total lifetime of the system does 
not change by more than 10% when going from SPS to RHIC energies. Finally, we have 
found in our model that chemical freeze-out of (anti-)baryons does not occur at the phase 
boundary and precedes the kinetic freeze-out of the system. 
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FIG. 1. Freeze-out time and transverse radius distribution d^N/^rxdrTdTfrdy) for pions (left 
column) and protons (right column). The top row shows the result for the pure hydro case up to 
hadronization with subsequent hadron resonance decays (but without hadronic reinteraction) . The 
bottom row shows the analogous calculation, but with full microscopic hadronic collision dynamics 
after the hadronization. The contour lines have identical binning within each column but differ 
between the two columns. 
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FIG. 2. Transverse freeze-out radius distributions d N/rTjdrxjdy for various hadron species. 
The distributions for vr, K, p, A and H are broad and similar to each other, whereas the 0,~ exhibits 
a narrower freeze-out distribution. 
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FIG. 3. Freeze-out time distributions d A^/drj,.dy of vr, p and i7^ for SPS and RHIC. Apart 
from the different integral values there is no significant difference between the RHIC and SPS 
distributions, i.e. the total lifetime of the reaction is comparable. 
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FIG. 4. Top: time evolution of on-shell hadron multiplicities (integrated over tt)- The dark 
grey shaded area shows the duration of the QGP phase whereas the light grey shaded area depicts 
the coexistance phase. Bottom: Hadron-hadron collision rates. 
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